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WRF 模式 在 天 山地 区 模拟 能 力 的 敏感 性 评估 
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摘 E: 为 探寻 天 山地 区 气候 模拟 WRE ( weather research and forecasting) 模式 中 最 优 参数 化 方案 组 合 , 针 对 云 微 物 
理 方 案 ( MIC) 、 积 云 对 流 方案 (CS) ,行星 边界 层 / 近 地 面 层 方案 (PBL/SLS) 、 陆 面 过 程 方案 (LSM ) 以 及 长 短波 辐射 


方案 (LSW ) 设 计 了 6 组 季节 尺度 物理 参数 化 方案 敏感 性 


试验 ,模拟 时 间 设 为 2014 -11 -28 一 2015 -12 -01。 利 用 


地 面 气象 站 观测 数据 和 GPM (global precipitation measurement) 卫星 降水 数据 (R=0.6) ,对 模式 模拟 的 日 最 高 .最低 


气温 及 降水 进行 验证 。 结 果 表 明 :WRF 模式 对 气温 的 模拟 效果 较 好 , 且 对 日 最 高 气温 (0.8 < R < 0.95) 的 模拟 优 于 


日 最 低 气 温 (0.62 < R «0. 88) ,日 最 低 气 温 对 陆 面 过 程 方案 和 云 微 物理 方案 敏感 。 模 式 对 降水 (R=0.6) 的 模拟 效 
果 一 般 , 各 方案 间 表 现 差异 大 。 云 微 物理 方案 WSM 6 - class 计算 量 较 小 且 能 较 好 地 模拟 降水 , 积 云 对 流 方案 Kain- 


Fristch 无 法 准确 处 理 天 山地 区 强 对 流 与 复杂 地 形 ,夏季 模拟 出 现 日 最 低 气 温 极 大 偏差 和 负 值 降水 。 得 出 的 最 优 组 
合 为 : WSM 6 - class ( MIC) , Betts-Miller-Janjie scheme ( CS ) , Mellor-Yamada-Janjic scheme/Monin-Obukhov ( Janjic 
Eta) , scheme( PBL/SLS) ,NOAH ( LSM) „Community Atmosphere Model ( LSW) 。 
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WRF 模式 (weather research and forecasting mod- 
el) "是 为 大 气 研究 和 业务 预报 应 用 而 设计 的 新 一 
代 中 尺度 数值 天 气 预 报 系统 ,完全 可 压 的 非 静 力 模 
式 , 被 广泛 应 用 于 几米 到 几 百 公里 尺度 的 各 类 气象 
气候 研究 中 。WRF 模式 提供 了 多 套 物理 参数 化 方 
案 ,不 同 尺 度 .地 理 位 置 以 及 应 用 目标 的 研究 ,所 需 
的 物理 参数 化 方案 组 合 各 不 相同 。 因 此 ,在 开展 区 
域 气候 模拟 之 前 ,需要 根据 特定 的 研究 区 域 和 研究 
目标 优选 模式 物理 参数 化 方案 组 合 。 

在 过 去 的 33 a(1979 一 2011 年 ) ,中 亚 地 区 的 气 
温 以 每 10a0.36 ~0.42 的 速率 在 上 升 中 ,天 山 是 
中 亚 地 区 的 “水 塔 " 中 ,气温 上 升 加 剧 了 天 山 雪 冰 的 
消融 ”” ,尽管 早期 会 带 来 充沛 的 雪 冰 融 水 ,但 在 其 
他 来 源 不 变 的 情况 下 ,冰川 体积 减 小 最 终 会 导致 其 
补给 径流 总 量 的 减少 ,从 而 影响 绿洲 区 人 民 的 生产 
生活 。 中 亚 是 “丝绸 之 路 经 济 带 ”的 核心 区 ,在 天 山 
开展 气息、 生态 、 水 文 等 研究 能 为 “丝绸 之 路 经 济 
带 ” 的 建设 提供 科学 支持 。 
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数据 是 科学 研究 的 灵魂 。 现 存 的 气候 数据 ,无 
论 是 全 球 气 候 模 式 general circulation model, 简称 
GCM) 数 据 或 是 各 类 再 分 析 数 据 , 水 平分 辩 率 都 在 
百 公 里 尺度 ,无 法 满足 高 精度 研究 需求 。 虽 有 少量 
高 分 辩 率 气候 数据 集 , 如 水 平分 辩 率 在 1 km 左右 
的 WorldClim 34; ( http ;// www. worldclim. org/) ,但 
其 缺少 垂直 高 度 场 数据 、 时 间 分 状 率 低 ( 月 尺度 )、 
提供 变量 较 少 [ 最 高 ( 低 ) 气 温 、 平 均 气 温 、 降 水 、 太 
阳 辐 射 、 风 速水 汽 压 ] 等 因素 , 仍 无 法 满足 高 时 间 
分 辨 率 多 维 数据 的 需求 。 利 用 区 域 气候 模式 (re- 
gional climate model, ,简称 RCM ) 对 现存 的 粗 分 辩 率 
气候 数据 , 即 GCM 数据 或 再 分 析 数 据 进行 动力 降 尺 
EE ,得 到 较 高 分 辩 率 的 气候 数据 集 是 常用 的 获取 高 
DRR 多维 气 候 数 据 集 的 方法 '" 。RCM 通过 内 部 
物理 过 程 作用 ,不 仅 可 以 保留 驱动 数据 的 大 尺度 特 
征 , 还 能 模拟 出 更 多 区 域 尺度 的 气候 特征 , 即 “ 增 
O, WRF 模式 作为 新 一 代 区 域 气候 模式 , 正 被 
越 来 越 多 的 研究 人 员 用 于 动力 降 尺 度 。 
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WRF 模式 物理 参数 化 方案 敏感 性 研究 在 多 地 
开展 , Argüeso 等 加 探究 了 不 同 WRF 参数 化 方案 对 
南西 班 牙 10 a 区 域 气候 模拟 的 敏感 性 发 现 ,日 最 低 
气温 对 不 同 的 行星 边界 层 方案 敏感 , 积 云 对 流 方案 
和 行星 边界 层 方案 对 降水 有 重要 影响 ,而 更 复杂 的 
云 微 物理 方案 则 没有 显著 影响 。1990 一 1995 年 
Mooney 等 中 开展 了 欧洲 区 域 气候 WRF 参数 化 方案 
的 敏感 性 评估 ,测试 了 2 组 长 波 辐射 方案 ,2 组 陆 面 
过 程 方案 2 组 云 微 物理 方案 以 及 2 组 行星 边界 层 
方案 ,发 现 2 m 气温 对 陆 面 过 程 方 案 最 敏感 、 云 微 物 
理 和 行星 边界 层 方案 次 之 ,并 且 ,WRF 不 能 很 好 地 
模拟 当地 降水 ,模式 结果 与 观测 相关 性 很 差 ,相关 系 
数 仅 为 0.1 ~0.3, Crétat °° FE WRF 物理 参数 方 
案 敏 感性 对 南非 区 域 气候 模拟 的 不 确定 性 进行 了 人 研 
究 ,模拟 时 间 设 为 1993 - 12 一 1994 - 02 南非 雨季 ， 
设计 了 由 3 个 积 云 对 流 方 案 .行星 边界 层 方案 以 及 
云 微 物 理 方案 组 成 的 27 组 试验 ,WRF 准确 刻画 了 
季节 性 大 尺度 降雨 格局 .降雨 梯度 以 及 气温 和 湿度 
的 垂直 廓 线 ,所 有 参数 方案 对 降雨 量 最 敏感 , 积 云 对 
流 方案 对 降雨 的 地 理 位 置 .强度 和 季节 指标 最 敏感 ， 
行星 边界 层 方案 和 云 微 物理 方案 次 之 。Flaounas 
等 加) 探究 了 WRF 模式 的 积 云 对 流 方案 和 行星 边界 
层 方案 敏感 性 对 人 们 理解 西非 季风 机 制 的 影响 ,发 
现行 星 边 界 层 方案 对 温度 、 湿 度 的 垂直 分 布 和 降雨 
量 有 和 较 大 影响 ,而 积 云 对 流 方案 对 一 些 动态 以 及 与 
降雨 相关 的 变量 影响 更 大 。Evans 等 ”在 澳洲 东海 
岸 针 对 一 系列 降雨 事件 进行 的 36 组 多 物理 参数 化 
方案 敏感 性 试验 显示 ,没有 一 组 组 合 方案 能 对 所 有 
事件 .变量 和 评价 指标 表现 最 佳 ,表明 不 同 的 气候 因 


70°E 75° E 80° E 


ChinaXiv 合 作 期 刊 


Wr 究 


子 对 不 同 的 物理 参数 方案 敏感 ,MYJ 行星 边界 层 方 
RF BMJ 积 云 对 流 方案 表现 稳定 ,不 建议 组 合 使 用 
YSU 行星 边界 层 方 案 、KF 积 云 对 流 方案 和 RRTMG 
长 短波 辐射 方案 ,同时 指出 ,降雨 强度 增加 时 物理 参 
数 化 方案 的 选择 变 得 尤为 重要 。Kala 等 "评估 了 
西 澳洲 西南 部 季节 尺度 WRF 对 驱动 数据 和 物理 方 
案 的 敏感 性 得 出 ,辐射 方案 的 选择 对 气温 和 降水 都 
有 较 大 影响 ,日 最 低 气 温 对 行星 边界 层 方案 十 分 敏 
感 , 积 云 对 流 方案 或 更 复杂 的 云 微 物理 方案 对 降水 
模拟 并 没有 显著 效果 。 

本 研究 以 天 山地 区 为 例 ,对 WRF 模式 物理 参 
数 化 方案 在 该 地 区 模拟 能 力 敏 感性 进行 了 评 佑 , 力 
图 找 出 一 套 适合 该 地 的 WRF 模式 物理 参数 化 方案 
最 优 组 合 ,为 该 地 区 气候 模拟 研究 提供 参考 。 


1 研究 区 概况 


天 山 是 亚洲 中 部 最 大 的 山系 , 西 起 乌兹别克 斯 
坦 的 克 孜 尔 库 姆 沙漠 以 东 , 经 塔吉克 斯 坦 哈萨克 斯 
坦 和 吉尔 吉 斯 斯 坦 进 入 我 国 新 疆 境 内 , 渐 失 于 哈密 
市 以 东 的 戈壁 中 (图 1)。 山 区 地 形 复杂 ,平均 海拔 
4 000 m 左右 ,最 高 的 托 木 尔 峰 达 7 435.3 m, KM 
横向 长 度 超 过 2 500 km ,南北 宽度 一 般 为 250 ~ 350 
km, 帕 米尔 以 北 的 天 山 山体 最 宽 , 可 达 800 km 以 
.EUU 。 由 于 巨大 的 高 度 和 有 利 的 地 理 位 置 以 及 独 
特 的 山地 走势 ,天 山 截留 了 大 量 西 风 环 流 所 携带 的 
大 西洋 水 汽 和 冬季 西北 风 所 携带 的 北冰洋 水 汽 , 夏 
季 南 亚 高 压 系统 发 达 , 天 山北 部 纬 向 环流 强大 ,形成 
天 山 夏 季 的 十 区 。 故 山区 积聚 了 丰富 的 降水 、 积 雪 、 
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图 1 天 山地 区 地 形 及 气象 观测 站 点 分 布 示 意 医 
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Fig.1 Topographical situation in the Tianshan Mountains and spatial distribution of the meteorological stations 
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冰川 等 形式 的 水 资源 ,中 亚 地 区 众多 河流 .湖泊 均 发 
源 于 此 ,哺育 了 我 国 最 长 的 内 陆 河 一 一 塔里木 河和 
水 量 最 大 的 内 陆 河 一 一 伊犁 河 ,是 西北 干旱 区 的 
"E" ,更 是 中 亚 地 区 的 “水塔 ”。 


2 试验 设计 与 评估 


2.1 模式 配置 

选用 WRF 模式 的 欧洲 中 期 天 气 预报 中 心 (Eu- 
ropean Centre for Medium-Range Weather Forecasts ) 
所 提供 的 再 分 析 数 据 ERA-Interim ^ zl Hy Bie BE; 
WRF 模式 为 大 气 物理 过 程 提 供 了 多 种 参数 化 方案 ， 
主要 包括 云 微 物理 方案 (cloud-microphysics, 简称 
MIC) 、 积 云 对 流 方案 (cumulus/convective scheme , 简 
Bk CS) .行星 边界 层 方 案 ( planetary boundary layer, 
简称 PBL)/ 近 地 面 层 方案 (surface layer scheme, 简 
TK SLS) , lil; i iE EE 77 3€ (land surface model, 简称 
LSM) 以 及 长 短波 辐射 方案 (longwave and shortwave 
radiation ,简称 LW and SW) ey 

选用 WRF V 3.9.1, 模 拟 区 域 ( 图 2) 采 用 双 层 
HAE (one-way nesting) ,粗细 网 格 水 平方 向 格 距 
分 别 为 40 km 和 8 km, 格 点 数 100 x 50,341 x 136, 
网 格 中 心 点 80" 下 ,42"N。 模 式 在 垂直 方向 上 设 有 
32 层 , 最 高 层 固 定 在 530 hPa。 以 ERA-Interim 数据 
集 ( 水 平分 辩 率 0. 75? x 0. 75°) 作为 模式 运转 的 初 
始 场 和 边界 场 ,边界 条 件 每 6 h 更 新 1 次 。 外 层 开 
启 网 格 松弛 ( grid nudging)"”, 系数 设 为 2 x 1077, 
用 于 松弛 模式 的 风 .气温 和 水 汽 ,使 其 忠实 于 驱动 条 
件 ,从 而 避免 在 长 时 间 模 拟 中 发 生 模 式 漂移 ( model 
drift) 。 内 层 关闭 网 格 松弛 ,体现 WRF 模式 动力 降 
尺度 的 “增值 ”效果 。 模 拟 区 内 ,各 问 边 界 上 设置 5 
个 格 点 的 边界 场 缓冲 区 (relaxation area ) ,使 缓冲 区 
内 模式 模拟 值 接近 大 尺度 强迫 因子 ,同时 与 模式 运 
行 得 出 的 中 尺度 数值 有 效 混合 。 开 启 海 洋 和 湖泊 ， 
如 巴尔 喀什 湖 (Balkhash lake ), 伊 塞 克 湖 ( Issyk 
Kul) 等 ,表面 温度 (sea surface temperature ) 实时 更 
新 。 采 用 WRF 模式 提供 的 USGS( United States geo- 
logical survey ) 地 形 资料 、MODIS ( moderate resolution 
imaging spectroradiometer ) 土地 利用 资料 、 AVHRR 
(advanced very high resolution radiometer ) 植被 覆盖 
度 资 料 以 及 FAO (food and agriculture organization of 
the United Nations ) 土壤 类 型 资料 。 

敏感 性 试验 模拟 时 间 段 为 2014 - 11 -28 一 
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图 2 WRF 模拟 区 域 
Fig.2 WRF domains 


2015 -12 - 01 ,模拟 时 间 前 3 d 将 作为 模式 模拟 的 


场 缓冲 区 的 模拟 区 域 ,将 被 用 于 模式 结果 验证 。 
2.2 试验 设计 

采用 Qiu 等 中 亚 地 区 动力 降 尺度 研究 所 选 
用 的 物理 参数 化 方案 作为 参考 组 (REFerence) , 即 : 
云 微 物理 方案 Thompson? 、 积 云 对 流 方 案 Betts- 
Miller-Janjic scheme!” (BMJ)、 陆 面 过 程 方案 No- 
ah?!) 行星 边界 层 方案 Mellor-Yamada-Janjic 
scheme?) ( MYJ ) 和 近 地 面 层 方案 Monin-Obukhov 
(Janjic Eta) scheme( MO) 长 短波 辐射 方案 Commu- 
nity Atmosphere Model?) (CAM)。 修 改 各 个 物理 参 
数 化 方案 进行 敏感 性 试验 ,试验 设计 详 见 表 1 。 

数值 模式 中 将 湿 物理 过 程 分 为 显 式 云 方案 ( 云 
微 物理 方案 ) 和 积 云 对 流 方 案 。WRF 中 云 微 物理 方 
案 主要 处 理 网 格 尺度 的 降水 ,试验 WSM6 所 采用 的 
WSM 6 - class 云 微 物理 方案 与 REF 所 使 用 的 
Thompson 方案 都 能 在 WRF 模拟 中 表达 出 水 的 6 种 
形态 , 即 云 水 和 冰 十 E KAM (cloud water and 
ice, rain, snow, vapour and graupel) ,这 是 模拟 天 山 
山区 强烈 对 流 及 复杂 降水 形式 所 必 不 可 少 的 。 
WRF 中 其 他 的 3 - class, 5 -class 云 微 物理 方案 都 


表 1 敏感 性 试验 方案 


Tab.1 Programs of the sensitivity experiments 


se, MI UAR MUR Rid gO 
ENR WIR 射 方案 流 方案 面 层 方案 方案 
REF NOAH CAM CAM BMJ MYJ/MO Thompson 
KF NOAH CAM CAM KF MYJ/MO Thompson 
RTG NOAH RRTMG RRTMG BMJ MYJ/MO Thompson 
RUC RUC CAM CAM BMJ MYJ/MO Thompson 
WSM6 NOAH CAM CAM BMJ MYJ/MO  WSM 6 - class 
YSU NOAH CAM CAM BMJ YSU/MMS MO Thompson 


TK 


无 法 准确 表达 。WSM 6 - class 方案 是 单 分 量 方案 ， 
事先 定义 好 水 成 物 函 数 供 模 式 计 算 使 用 ,而 Thomp- 
son 方案 是 双 分 量 方案 ,预报 变量 包含 以 上 6 种 水 
成 物 的 比 水 含量 和 冰晶 数 浓度 ,更 接近 真实 情况 ,也 
EAS AROY 。 高 分 辩 率 天 气 过 程 的 模拟 往往 需要 较 
复杂 的 云 微 物理 方案 ,区 域 气候 模拟 则 可 放宽 。 
因此 ,有 必要 测试 对 比 复杂 与 略 不 复杂 但 计算 量 小 
的 方案 模拟 结果 ,在 模式 计算 量 与 模式 表现 间 寻 求 
平衡 。 

KF 与 REF 试验 相同 ,除了 将 REF 中 的 积 云 对 
流 方案 BMJ 改 为 Kain-Fristch" (KF), WRF rp fi 
云 对 流 方案 用 于 处 理 次 网 格 尺度 深 对 流 和 浅 对 流 积 
云 产 生 的 大 尺度 效应 ,并 描述 网 格 无 法 分 辨 的 云 内 
上 升 气流 下降 气 流 以 及 云 外 补偿 气流 ,只 有 在 单个 

让 云 中 被 触发 时 才 起 作用 , 即 调整 垂直 加 热 和 湿度 

廓 线 并 产生 对 流 降 水 。BMJ 方案 由 Betts-Miller 77 
案 改 进 获 得 ,使 深 对 流 特 征 廓 线 及 松弛 时 间 可 随 积 
云 效 率 发 生变 化 , 浅 对 流水 汽 特 征 廊 线 中 炉 的 变化 
较 小 且 为 非 负 值 。KF 方案 利用 一 个 简单 的 云 模式 
表达 气流 抬升 .下 沉 和 卷 人 、 卷 出 ,是 WRF-ARM 
(advanced research WRF ) 推荐 使 用 的 气候 模拟 积 云 
对 流 参数 化 方案 "" 。 

行星 边界 层 方案 Yonsei University ( YSU ) 
schem 7? gi HERI YSU 中 替代 REF 中 的 MYJ 
方案 ,由 于 WRF 中 行星 边界 层 方案 与 近 地 面 层 方 
案 一 一 对 应 , 故 近 地 面 层 方案 也 由 MO 改 为 MMS 
Monin-Obukhov scheme( MM5 MO), KAHH JZP 
Hm ER KAA E .水汽 及 动量 进行 垂直 输送 TEL 
合 , 故 行星 边界 层 能 确定 混合 层 和 稳定 层 内 的 通 量 
FRAG ,提供 整个 空气 柱 内 的 大 气温 度 .水汽 及 动量 通 
量 倾向 。MYJ 方案 由 2 阶 注 流 闭合 方案 简化 而 来 ， 
采用 TKE (turbulent kinetic energy ) 闭合 技术 ,在 边 
界 层 内 部 预报 满 流动 能 ,对 所 有 的 灌流 阶 量 进行 诊 
断 , 最 终 闭合 边界 层 内 动量 方程 ”) 。YSU 方案 由 
MRF( Medium Range Forecast) 行 星 边 界 层 方案 改进 
而 来 ,采用 下 闭合 技术 ,增加 了 对 边界 层 顶 卷 层 的 
显 式 处 理 , 解 决 了 原 MRF 方案 中 过 度 混合 的 缺点 ， 
也 是 WRF-ARM 推荐 使 用 的 气候 模拟 方案 ” 。 

RUC 试验 主要 用 于 测试 陆 面 过 程 方案 对 模式 
模拟 结果 的 敏感 性 , 陆 面 过 程 方 案 能 显著 有 影响 陆 一 
气 间 动 量 、 热 量 和 水 汽 交 换 。REF 试验 所 采用 的 
Noah 方案 是 WRF 区 域 气候 研究 中 最 常用 的 陆 面 过 
程 方案 ,有 4 个 土壤 层 :0 ~10 em,10 ~40 cm 40 ~ 


100 cm 、100 ~ 200 em, 本 次 试验 中 的 RUC FR 
与 Noah 同等 复杂 ,设置 6 个 土壤 层 :0 cm、5 cm 、20 
em ,40 cm、160 cm 、300 cm ,但 尚未 被 充分 评 佑 使用。 
试验 RTG 与 REF 的 不 同 之 处 是 采用 了 RRT- 
MG ”作为 长 短波 辐射 方案 ,RRTMG 是 短波 辐射 方 
案 RRTM (rapid radiative transfer model) 应 用 在 GCM 
模拟 时 的 修正 版 。 长 短波 辐射 的 精确 刻画 对 低层 气 
温 模 拟 十 分 关键 ,CAM 与 RRTMG 长 短波 辐射 方案 
相 比 ,WRF 模拟 中 常见 的 RRTM/Dudhia 长 短波 辐 
射 方案 组 合 更 为 复杂 ,CAM 与 RRTMG 都 考虑 到 了 
云层 重奏 的 情形 ,以 判断 网 格 云 量 ,而 后 者 只 是 以 二 
进 制 的 方式 简单 测算 是 否 有 云 。CAM 与 RRTMG 
还 把 痕 量 气体 气 深 胶 、 臭 气 和 二 氧化 碳 效应 ,甚至 
反射 的 短波 辐射 通 量 都 考虑 在 模式 计算 中 ,但 CAM 
使 用 Delta-Eddington 估算 短波 辐射 的 吸收 与 散 
gr?» , RRTMG 则 是 用 Correlated-k 辐射 传输 算 
法 ,试验 RTG 和 欲 测 试 RRTMG 长 短波 辐射 方案 能 
否 提 高 模式 在 天 山 模拟 的 精度 。 
2.3 ”验证 数据 

利用 地 面 气象 观测 站 的 观测 资料 ,对 模式 模拟 
的 日 最 高 .最低 气温 及 降水 进行 验证 。 笔 者 所 使 用 
的 地 面 气象 资料 来 源 于 中 国 气象 数据 网 (http ://da- 
ta. cma. cn/ ) 以 及 GHCN (global historical climatology 
network) daily observation data map ( https ://gis. nc- 
de. noaa. gov/maps/ncei/cdo/ daily ) °°? 所 提供 的 逐 
日 最 高 气温 .最低 气 温 数据 和 24 h 降水 量 。 研 究 中 
对 所 有 站 点 进行 质量 控制 ,对 缺 测 数据 大 于 20% 的 
站 点 直接 剔除 不 予 考虑 ,最 终 在 研究 区 范围 内 选 出 
131 个 气温 观测 站 .53 个 降水 观测 站 (图 1)。 

GPM ( global precipitation measurement ) 核心 卫 
Fe 2014 4E 2 月 27 日 发 射 , 是 新 一 代 降 水 观测 卫 
星 ,分 辨 率 达 0. 1°, 首次 采用 双 频 雷达 观测 系统 
(DPR) ,结合 主动 雷达 ,从 不 同 角度 提供 云层 降水 
粒子 的 物理 信息 ,如 雨滴 的 形状 、 强 度 以 及 对 流 过 
程 ,提高 了 对 轻 量 级 降水 的 捕捉 能 力 ,能 更 精准 地 刻 
画 降 水 的 时 空 分 布 ” 。 金 晓 龙 等 在 天 山地 区 对 
比 GPM 卫星 和 另外 两 套 应 用 较为 广泛 的 卫星 降水 
产品 TRMM 与 CMORPH 发 现 ,CPM 卫星 对 天 
山 山区 降水 捕 提 表现 最 佳 (R=0.6)。 本 研究 利用 


GPM 卫星 降水 数据 ,对 内 层 模 拟 的 降水 结果 进行 
验证 。 
2.4 评估 方法 


人 研究 中 对 各 个 试验 在 冬 ( December-January-Feb- 
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ruary , DJF) 、 春 (March-April-May,MAM) , € ( June-Ju- 
ly-August, JJA ) 、 秋 ( September-October-November , 
SON)4 个 季节 的 表现 进行 评估 ,评估 要 素 有 季节 平 
均 日 最 高 气温 .日 最 低 气 温 以 及 月 累积 降水 。 首 先 ， 
将 日 值 气象 站 点 数据 转换 到 季节 斥 度 。 研 究 选 用 最 
近邻 法 获得 气象 站 点 的 模拟 值 ,从 而 减少 其 他 插值 
MERAT fü oe ng veas 7 07 。 利 用 双 线 性 搬 值 
算法 ,对 天 山地 区 0. 1? 分 辩 率 月 值 GPM 卫星 降水 
数据 重 采样 至 8 km 水 平分 辨 率 ,同时 也 转换 到 季节 
尺度 。 根 据 实际 地 形 资 料 对 模拟 得 到 的 气温 数据 进 
行 地 形 效 应 校正 ,由 于 使 用 模式 自身 的 递减 率 与 
标准 递减 率 校 正 时 影响 不 大 ”9 , 故 研究 中 取 标 准 大 
气温 度 直 减 率 系数 0.65 名， (100m) ”作为 地 形 校 
正 系 数 。 

Taylor 图 是 从 两 个 场 的 空间 相关 系数 (R ) .中 
心 化 均 方 根 误差 (8') .模拟 场 的 标准 差 (o, ) 和 观测 
场 的 标准 差 (o,) 这 4 个 统计 量 来 比较 模式 模拟 结 
果 与 观测 场 的 吻合 程度 。 横 纵 坐 标 代 表 模 拟 场 
与 观测 场 的 标准 差 之 比 , 连 接 横 纵 坐 标 尾 端 的 圆 弧 
代表 模拟 场 与 观测 场 的 空间 相关 系数 , 横 轴 上 以 
REF 点 为 圆心 的 半圆 弧 是 中 心 化 均 方 根 误差 。 简 
TEL, Taylor 图 结果 越 靠 近 标 准 差 之 比 为 1 的 圆 弧 
线 (模拟 场 和 观测 场 的 标准 差 越 相近 , 即 空 间 变 率 
越 接 近 ) 、 越 接近 横 轴 (中 心 化 均 方 根 误 差 越 小 、 模 
拟 场 和 观测 场 的 空间 相关 系数 越 大 ) ,模拟 效果 越 
Xf, Taylor 图 已 用 于 气候 模式 的 性 能 评 
4d 2 795999 ,本 次 研究 利用 Taylor 图 对 模式 模拟 
的 多 组 试验 结果 进行 评价 。 但 Taylor 图 评估 的 是 模 
式 对 距 平 场 的 模拟 能 力 ,不 能 很 好 地 将 模式 与 观测 
的 平均 值 差异 反映 出 来 。 因 此 ,本 研究 还 将 利用 偏 
差 ,包括 绝 对 偏差 (absolute bias ,简称 AB) 和 百分比 
{iE (percentage bias, 简称 PB ) 来 对 模式 结果 进行 
评估 。 评 估 指 标的 计算 公式 如 下 : 


N = n 
R = (M, - M) (0; -O)/o,, 9, 


Pe 12 M n 2 
E = [XL -M)(0, -0)] 


N — — 
X(M,-0,)-M-0 


m (M, = 0;) 
25.0, 


式 中 :假定 待 评估 要 素 有 IN 个 模拟 值 (只 ) 和 观测 值 


PB = 


3 结果 分 析 


3.1 气温 

图 3 和 图 4 分 别 是 日 最 高 气温 和 日 最 低 气 温 
Taylor 图 ,所 有 试验 空间 相关 系数 都 通过 了 99. 9% 
的 置信 检验 ,偏差 结果 见 表 2 和 表 3。 各 个 试验 模 
拟 的 日 最 高 气温 空间 相关 系数 在 0. 8 ~ 0.95, RE 
4 RTG RUC 外 ,中 心 化 均 方 根 误差 均 控制 在 第 三 
级 ( 即 图 中 以 REF 点 为 圆心 向 外 的 第 3 个 半圆 弧 ) 
以 内 ,标准 差 比率 仅 有 秋季 YSU、RTG 及 夏季 位 于 
1+0.25 以 外 。 相 对 而 言 ,模式 对 日 最 低 气 温 的 模 
拟 效果 不 如 日 最 高 气温 。 日 最 低 气 温 Taylor K| w 
示 ,各 个 试验 标准 差 比率 接近 1( +0.25 左右 ) , ME 
有 试验 YSU 、RTG 夏季 模拟 除外 ,空间 相关 系数 在 
0.62 ~0. 88 ,中 心 化 均 方 根 误差 相 比 日 最 高 气温 模 
拟 普 遍 偏 大 。 

从 Taylor 图 上 看 ,REF 试验 对 日 最 高 气温 的 模 
拟 全 年 表现 较 好 、 较 稳定 ,秋季 最 佳 ,系统 性 低 佑 观 
测 场 均值 ,最 大 低估 出 现在 冬季 和 春季 ,偏差 达 
-3.06 CFU -3.33 C( - 36. 26% Fil — 12. 27%). 
REF 试验 对 日 最 低 气 温 的 模拟 全 年 也 较 稳定 ,但 除 
夏季 外 都 高 估 了 观测 场 均 值 ,最 大 绝对 偏差 出 现在 
冬季 ,高 估 1.75 "C ( -10.3% ) ,最 大 百分比 偏差 在 
秋季 , 达 -50.2% (1.3%)。 

日 最 高 气温 模拟 上 ,Taylor 图 中 各 个 试验 的 模 
拟 结果 均 较 好 上 且 相似 , 唯 夏季 RTG RUC ,冬季 YSU、 
RTG 表现 相对 较 差 。 偏 差 评 佑 中 ,试验 YSU EF, 
夏季 绝对 偏差 较 大 ,为 -5.89 CH -8.78 %C ,试验 
KF REF 在 冬季 百分比 偏差 较 大 ,为 -37.16% 和 
-36.26% 。 但 在 日 最 低 气 温 模拟 上 , Taylor 图 中 唯 
有 冬季 各 个 试验 的 表现 较 接 近 , 其 余 3 个 季节 里 ， 
REF , WSM6 , RUC 效果 较 好 ,剩余 试验 表现 较 差 。 
值得 一 提 的 是 ,试验 RUC 仅 在 夏季 略 逊 于 试验 
REF ,其 余 季 节 均 表现 最 佳 。 偏 差 评估 中 ,试验 KF 
(H: -890% ) RTG( Æ: -104.66%)、RUC( 秋 : 
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表 2 日 最 高 气温 季节 绝对 偏差 和 百分比 偏差 
Tab.2 Seasonal absolute and percentage biases 


of daily maximum temperature 


试验 REF KF RTG RUC WSM6 YSU 
DJF -3.06 -3.14 0.82 -1.38 -1.94 一 0.87 
-36.26% -37.16% 9.68% -16.39% -22.97% -10.32% 
MAM -3.33 -2.55 -2.07 -1.17 -3.09 -5.89 
-12.27% -9.39% -7.65% -4.30% -11.41% -21.70% 
JJA -0.69 =2.75 -4.46 2:29 -1.99 -8.78 
-1.90% -7.56% -12.25% 6.29% -5.46% -24.12% 
SON -1.64 -2.52 -1.67 1.36 -1.42 -4.15 
-6.83% -10.47% -6.93% 5.65% -5.91% -17.25% 


表 3 日 最 低 气温 季节 绝对 偏差 和 百分比 偏差 
Tab.3 Seasonal absolute and percentage biases 
of daily minimum temperature 
试验 REF KF RTG RUC WSM6 YSU 
DJF 1.75 3.43 6.62 3.39 3.66 4.38 
—10.3096 -20.14% -38.86% -19.90% -21.52% -25.70% 


MAM -0.79 -0.28 2.62 1.69 - 0.56 0.65 
31.36% 11.15% —104.6696 -67.59% 22.52% -26.00% 

JJA | 0.84 -92.82  -0.65 3.14 0.64 -3. 776 
8.04% -890.68% -6.28% 30.1296 6.09% -36.13% 

SON 1.3 0.46 2.06 2.75 0.67 -0.57 


-50.20% -17.62% -79.84% — 106.35% -26.01% 22.15% 


-106.3596 ) 都 出 现 了 超 100% 的 负 偏 差 , 尤 其 KF 
试验 在 模拟 过 程 中 极 有 可 能 存在 积分 溢出 的 情况 ， 
试验 WSM6 除 冬 季 偏 差 较 大 3.66( -21.52% ) ,其 
余 季 节 均 优 于 REF, 
3.2 降水 
图 5 是 月 累积 降水 气象 站 点 和 GPM 卫星 验证 
Taylor 图 ,偏差 评估 见 表 4 ~5。 除 试验 YSU( 春 、 夏 、 
BK) .KF( 夏 , 负 相 关 ) 、RTG( 夏 ) 外 ,大 部 分 试验 都 
通过 了 99.9% 空间 相关 系数 置信 检验 。 模 式 模拟 
的 降水 空间 相关 系数 除 冬 季 站 点 验证 (KF、YSU f 
差 ,REF 接近 0.9,WSM6 RTG „RUC 都 超过 了 0.9) 
外 都 普遍 偏 低 ,为 0.6 左右 或 更 差 ,但 与 Mooney 
等 ”在 欧洲 模拟 的 降水 效果 (0.1 <R <0.3) 相 比较 
好 。 标 准 差 比率 大 部 分 试验 都 在 1 +0.25 外 ,试验 
间 差 异 大 (秋季 差异 最 大 ) 。 中 心 化 均 方 根 误差 除 
冬季 表现 较 好 的 试验 (RUC .RTG , WSM6 , REF) 基本 
都 在 第 四 级 或 第 四 级 以 外 。 

试验 REF 对 降水 的 模拟 冬季 较 好 ,秋季 最 差 ， 
对 观测 场 空间 变 率 刻画 是 所 有 试验 中 最 好 的 。 系 统 


性 低估 隆 水 均值 。 站 点 验证 冬季 偏差 最 小 ,为 
-0.22( -3.88%); 最 大 偏差 出 现在 夏季 ,为 
-10.63( -36.21% ) ,卫星 验证 4 个 季节 偏差 都 在 
28% 左右 。 

站 点 验证 中 ,冬季 ,试验 WSM6,RTG,RUC 与 
REF 的 表现 接近 或 稍 好 。 春 季 ,RUC 5j WSM6 标准 
差 比率 劣 于 REF ,但 中 心 均 方 根 误差 与 空间 相关 系 
数 均 优 于 REF, HÆ REF 与 WSM6 的 表现 相当 ， 
RUC 略 差 。 秋 季 各 个 试验 间 差 异 大 ,REF 标准 差 比 
率 接近 于 1 ,但 中 心 均 方 根 误差 过 大 (第 五 级 ) .空间 
相关 系数 较 小 (0.3 左右 ) ,而 WSM6 以 接近 0. 85 的 
空间 相关 系数 以 及 相对 较 小 的 中 心 均 方 根 误差 (第 
三 级 ) ,虽然 标准 差 比 率 较 差 , 但 仍 可 判断 为 秋季 模 
拟 效果 最 好 的 试验 ,在 其 他 试验 普遍 出 现 失误 时 ， 
WSM6 表现 较为 稳定 。KF 试验 模拟 的 夏季 降水 出 
现 负 值 ,与 观测 场 存在 不 显著 负 相 关 ( -0. 067 ) ,这 
与 气温 (日 最 低 气 温 , 夏 ) 模拟 时 判断 该 试验 可 能 出 
现 积分 溢出 互 为 印证 ,其 他 季节 的 模拟 效果 也 较 差 。 
试验 YSU 在 多 个 季节 出 现 不 显著 相关 ,模拟 效果 较 
差 ,试验 RTG 在 冬季 表现 较 好 ,其 余 季 节 表 现 欠 


表 4 站 点 验证 的 月 累积 降水 季节 绝对 偏差 和 百分比 偏差 
Tab.4 Seasonal absolute and percentage biases of monthly 


accumulated precipitation validated at different stations 


试验 REF KF RTG RUC WSM6 YSU 
DJF -0.22 -4.43 -2.29 -0.3 -1.39 -4.24 
-3.88% -77.35% -39.91% -5.15% -24.20% -74.03% 
MAM -3.86 -13.51 -14.94 -6.21 -9.44 -16.69 
-22.76% -79.55% -87.97% -36.56% -55.57% -98.26% 
JJA -10.63 -27.69 -26.67 -11.06 -8.27 -28.50 
-36.21% -94.29% -90.82% -37.65% -28.16% -97.06% 
SON -3.31 -13.3 -0.18 3.52 -7.49 -11.33 
-20.84% -83.86% -1.15% 22.20% -47.25% -71.42% 


表 5 卫星 验证 的 月 累积 降水 季节 绝对 偏差 和 百分比 偏差 
Tab.5 Seasonal absolute and percentage biases of 


monthly accumulated precipitation validated by GPM 


试验 REF KF RTG RUC WSM6 YSU 

DJF -6.09 -16.25 -2.71 0.06 -2.03  -13.59 
-31.37% -83.77% -13.99% 0.3296 -10.45% -70.06% 

MAM -8.46 -21.36 -17.16 -7.69 -12.34 -26.41 


-30.33% -76.57% -61.52% -27.58% -44.24% -94.65% 
JJA -5.40 -26.21 -13.14 -2.75 -3.83 -17.23 
-25.27% - 122. 60% -61.45% -12.85% -17.92% -80.57% 
SON -8.03 -17.18 1.73 4.78 -10.31 -11.31 
-32.20% -68.91% 6.93% 19.16% -41.36% -45.38% 
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图 5 站 点 和 卫星 验证 季节 月 累积 降水 泰勒 图 


Fig.5 Taylor diagrams of monthly accumulated precipitation validated at different stations and by GPM 


佳 。 偏 差 评估 中 ,所 有 试验 都 低 佑 了 研究 区 观测 场 
的 平均 降水 ,4 个 季节 里 REF, WSMG 在 所 有 试验 中 
表现 较 好 且 较 稳定 ,REF 对 观测 场 均值 的 刻画 仅 在 
REMF WSM6 秋季 差 于 RTG, 

GPM 卫星 验证 与 站 点 验证 结果 分 布 较 相似 , 试 
验 空间 相关 系数 在 0.6 左右 ,标准 差 比 率 和 中 心 化 
均 方 根 误差 相 比 ,站 点 验证 结果 稍 差 且 试验 间 差 异 
Ko ZE REF 试验 表现 优 于 其 他 试验 ,春季 试验 
WSM6 及 RTG 的 效果 较 好 ,夏季 仍然 是 REF 试验 表 
现 较 好 ;KF 试验 与 站 点 验证 类 似 ,也 出 现 了 负 相 关 ， 
秋季 试验 WSM6 表现 最 好 ;RTG 试验 的 中 心 化 均 方 


其 余 季 节 一 般 。 
4 结论 与 讨论 


本 文 在 中 亚 动 力 降 尺 度 研 究 基础 上 ,设计 了 6 
组 参数 化 方案 敏感 性 试验 ,以 探寻 天 山地 区 气候 模 
W WRF 模式 中 最 优 参 数 化 方案 组 合 ,同时 为 天 山 
地 区 气候 模拟 研究 提供 参考 。 利 用 气象 站 点 观测 数 
据 、GPM 卫星 降水 数据 ,对 模式 模拟 的 日 最 高 .最 低 
气温 及 降水 进行 验证 ,主要 结论 如 下 : 

(1) 模式 对 气温 的 模拟 效果 较 好 , 且 对 发 生 于 
白天 的 日 最 高 气温 刻画 优 于 晚上 的 日 最 低 气 温 。 日 


根 误差 特别 大 ,甚至 在 图 幅 外 。 偏 差 评估 中 ,与 站 点 
验证 普遍 低估 观测 场 均值 不 同 ,有 试验 在 秋冬 季 高 
估 降 水 ,这 与 卫星 对 秋冬 季 轻 量 级 降水 监测 困难 相 
KOM, pU; REF, WSM6,RUC 的 效果 较 好 , 且 试 验 
REF 在 4 个 季节 表现 平稳 , WSM6 在 冬夏 季 优 于 
REF ,而 在 春秋 季 差 于 REF, RUC 在 冬季 表现 最 佳 ， 


最 高 气温 模拟 的 各 试验 (实际 是 各 物理 方案 间 ) 未 
见 明显 差异 ,日 最 低 气温 对 陆 面 过 程 方 案 及 云 微 物 
理 方 案 敏 感 。 积 云 对 流 参 数 化 方案 KF 对 日 最 高 气 
温 的 模拟 较 好 ,但 在 日 最 低 气 温 模拟 ( 夏季) 中 出 现 
极 大 负 偏 差 ,存在 积分 溢出 异常 。 

(2) 模式 对 降水 的 模拟 效果 一 般 , 大 部 分 试验 
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都 低估 了 天 山地 区 的 平均 降水 量 。 与 气温 模拟 不 
同 ,降水 模拟 结果 间 差 异 较 大 ,但 没有 表现 特别 好 的 
试验 。 参 考 参 数 化 方案 组 合 (REF) 对 降水 的 模拟 尚 
可 , 云 微 物理 方案 WSM 6 - class 与 观测 场 均值 的 模 
拟 相 比 ,REF 稍 差 ,但 总 体 表现 最 佳 。KF 试验 在 模 
拟 夏季 降水 时 出 现 负 值 , 且 与 观测 场 有 一 个 不 显著 
的 负 相 关 ,表明 积 云 对 流 方案 KF 在 模拟 天 山地 区 
时 不 稳定 ,出 现 异 常 。 

(3) 没有 试验 在 各 项 评估 指标 中 均 表 现 最 优 ， 
这 与 Evans 42°”) 的 研究 结果 一 致 ,但 所 有 试验 中 
WSM6 综合 表现 较 好 。 同 时 , 单 分 量 云 微 物理 方案 
WSM 6 - class 比 双 分 量 Thompson 方案 在 计算 量 上 
有 所 减 小 ,这 能 为 高 计算 量 区 域 气候 模拟 减负 ,可 选 
取 试 验 WSM6 所 代表 的 参数 化 方案 组 合 , 即 
WSM 6 - class ( MIC) , BMJ ( CS) , MYJ/MO ( PBL/ 
SLS) .NOAH ( LSM) ,CAM ( LSW) 进行 天 山地 区 气候 


模拟 研究 。 

本 研究 主要 是 为 了 寻求 区 域 气候 模拟 最 优 参数 
化 方案 组 合 ,对 格 点 代表 性 误差 中 简单 处 理 ,后 期 
气候 模拟 研究 中 需 考虑 如 何 将 格 点 代表 性 误差 最 小 
化 。 降 水 观测 站 点 主要 集中 在 中 国 新 疆 境 内 的 东 天 
LL, ARA AA m, CRIES S| A CPM 卫星 降 
水 数据 验证 ,卫星 验证 结果 与 站 点 验证 差异 不 大 。 
积 云 对 流 方案 KF 被 广泛 应 用 在 WRF 气候 模拟 
中 -1 ,但 在 本 次 研究 中 发 生 积分 溢出 异常 , 即 夏 
季 模 拟 出 现 日 最 低 气 温 极 大 偏差 和 负 值 降水 ,该 方 
案 可 能 难以 处 理 天 山地 区 强 对 流 与 复杂 地 形 ,尤其 
在 环境 温度 较 高 ,对流 强烈 的 夏季 ,天 山地 区 区 域 气 
候 模 拟 不 宜 选 用 该 方案 ,考虑 调整 该 方案 内 部 的 垂 
直 温 湿 廊 线 , 如 设置 极 值 ,以 防 出 现 积分 溢出 、 无 法 
收敛 的 情况 。 同 时 也 表明 ,进行 WRF 区 域 气候 模 
拟 研究 时 ,首先 进行 物理 参数 化 方案 敏感 性 评估 是 
非常 必要 的 。 
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Abstract; In this study ,a sensitivity analysis of the weather research and forecasting model ( WRF) to different 
physical options in the Tianshan Mountains was carried out to seek the best performance combination. Yearly simu- 
lations at seasonal scale from 28th November ,2014 to 1st December ,2015 were carried out. The station observations 
and GPM (Global Precipitation Measurement) satellite precipitation data ( Rz20. 6) were used to validate the simu- 
lated daily extreme temperatures and precipitation. This work showed that temperature was modeled by WRF well, 
the simulated results of daily maximum temperature ( T2max,0. 8 < R «0.95) were better than that of daily mini- 
mum temperature ( T2min,0. 62 < R «0. 88) , and T2min was sensitive to land surface model and cloud-microphysi- 
cal model. The simulated results of precipitation by WRF were not so ideal with low correlation coefficients about 
0.6 (R«0. 6) , and they were different from different schemes. The cloud-microphysical scheme WSM 6 - class 
could be used to simulate precipitation well. Kain-Fristch cumulus/convective scheme could not be used to describe 
the strong convection and complex terrain in the Tianshan Mountains as we got the very large T2min biases and neg- 
ative values of precipitation in summer. The best performance combination was WSM 6 - class ( Cloud-microphys- 
ics) , Betts-Miller-Janjic ( Cumulus/Convective scheme ) , Mellor-Yamada-Janjic scheme/Monin-Obukhov ( Janjic 
Eta) scheme (Planetary boundary layer/Surface layer scheme) ,NOAH (Land surface model) ,and Community At- 
mosphere Model ( Longwave and shortwave radiation ) . 
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